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Abstract 

When a two-dimensional curved surface is conceived as a limiting case of a curved shell of equal 
thickness d, where the limit d — > is then taken, the well-known geometric potential is induced 
by the kinetic energy operator, in fact by the second order partial derivatives. Applying this 
confining procedure to the momentum operator, in fact to the first order partial derivatives, we 
find the so-called geometric momentum instead. This momentum is compatible with the Dirac's 
canonical quantization theory on system with second-class constraints. The distribution amplitudes 
of the geometric momentum on the spherical harmonics are analytically determined, and they are 
experimentally testable for rotational states of spherical molecules such as Cqq. 

PACS numbers: 02.40.-k, 03.65.-w, 68.65.-k, 73.22.Dj 



1 



Introduction The affirmative experimental evidence in 2010 i of the geometric potential 
firstly explored in 1971 2 and fundamentally finished in 1981 3 and with correct inclusion 
of electromagnetic field in 2008 1 etc. 5 is a groundbreaking advance of quantum mechanics 
applied for curved nanostructures, starting from the three dimensional (3D) bulk system and 
then reducing it to a 2D surface one. 2-5 This success echoes a historical footnote in Dirac's 
Principle on the canonical quantization assumption that "is found in practice successful only 
when applied with the dynamic coordinates and momenta referring to a Cartesian system of 
axes and not to more general curvilinear coordinates.'"' However, the classic work by Jensen, 
Koppe 2 and da Costa 5 hides an important physical and mathematical message when dealing 
with derivatives on a 2D curved surface S: There is a noninterchangeability of order of taking 
two limits. Explicitly, when the 2D curved surface is conceived as a limiting case of a curved 
shell of equal thickness d, where the limit d — > is then taken, great discrepancies present 
as firstly taking limit d — > then defining the derivatives on the surface, and as firstly 
defining derivatives in bulk then letting d — > 0. The second order is named as the confining 
procedure for studying motion on 2D surface embedded in 3D. 1 For the former order, the 
quantum kinetic energy operator is hypothesized to be proportional to Laplace-Beltrami 
operator Alb on the surface:' 

whereas for the latter order, we have following form of the kinetic energy operator T 
instead: 2-4 

r = < 2 > 

where M is the mean curvature and K is the gaussian curvature, and the excess term is 
called geometric potential V gp , [ ' 1 

V gp = -^(M*-K). (3) 

The experimental verification of this potential implies that the original Laplace-Beltrami 
operator Alb on the 2D surface may not be enough unless a term (M 2 — K) is included, 1 ' 8 

A LB ^ A LB + (M 2 -K). (4) 



For avoiding confusion, we adhere to the convention in mathematics where the Laplace 
operator V 2 acting on a function is defined by the divergence of the gradient of the function 
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in flat space: V 2 = V • V, while the Laplace-Beltrami operator Alb is a generalization of 
the Laplace operator on surface under consideration. 

The noninterchangeability of calculus order must be fundamentally associated with the 
gradient V, or the momentum operator p = — iftV, which has not explored before. On the 
other hand, we must mention an entirely independent development on the quantization of 
the momentum on 2D surface embedded in 3D flat space,' 1 and the momentum is found to 
assume the following form, 

p^-i^ + Mn), (5) 

where we use the tensor covariant and contravariant components and the Einstein summation 
convention, and 

r(q\q 2 )=(x(q\q 2 ),y(q\q 2 ),z(q\q 2 )) (6) 

is the position vector on the surface S parametrized by (q 1 ^ 2 ) denoted by (f and q v with 
lowercase greek letters /i, v taking values 1,2, and r M = g^ v y: v = g^ u d u r = g^dr/q 1 ' with 
5V = ' ^ r being the metric tensor. At this point r, n = (n x ,n y ,n z ) is the normal and 
Mn symbolizes the mean curvature vector field, a geometric invariant. 9 

The first aim of the present study is to show that the application of the same confining 
procedure pioneered by Jensen, Koppe 2 and da Costa' to operator p = —iftV that holds 
true in bulk automatically results in p = —ik^d^ + Mn) (5) on the surface. So, analogue 
to the name geometric potential we can call p = —ih{v^d^ + Mn) (5) geometric momentum 
(GM). Because r(g 1 ,g 2 ) (6) in mathematics offers the so-called standard parametrization of 
the 2D surface, the corresponding GM (5) should offer proper description of the momentum. 
If simply denoting the gradient operator r^d^ 10 on the surface by V//, Eq. (5) implies 
following correspondence, 

V// V// + Mn. (7) 

On a surface, is there a component Mn normal to it? This result (7) is somewhat contrary 
to what physical intuition or common sense would indicate. But it is the case as examined 
in 3D flat space. 

GEOMETRIC MOMENTUM AS A CONSEQUENCE OF CONFINING PROCEDURE To prove 

the GM (5), we utilize exactly the same manner how to derive the geometric potential. 2-4 
For ease of the comparison, we use similar set of symbols as Ferrari and Cuoghi who recently 
build up a theoretical framework with geometric potential when the electromagnetic field is 
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present. 4 The lowercase Latin letters k stand for the 3D indices and assume the values 
1,2,3, e.g., (xi,pj) for the position and momentum. Position specified by (g 1 ,^ 2 ,^ 3 ) can 
be understood as description of the position in the curvilinear coordinates parameterizing a 
manifold. The original 2D surface r(g 1 ,g 2 ) is considered as a more realistic 3D shell whose 
equal thickness d is negligible in comparison with the dimension of the whole system. The 
position R within the shell in the vicinity of the surface S can be parametrized as with 
< g 3 < d, 

R(q\ q\ q 3 ) = r(q\ q 2 ) + q 3 n(q\ q 2 ). (8) 

The gradient operator V in 3D flat space, expressed in the curvilinear coordinates, takes 
following form, 10 

V = r^ + n9 g 3. (9) 
The relation between the 3D metric tensor and the 2D one is given by, 3,1 

Gij = g^u + [ag + [ag) T ] ^ q 3 + {aga T ) yw (q 3 ) 2 , 

G M3 = G 3fl = 0, G 33 = 1, (10) 

where a^ u is the Weingarten curvature matrix for the surface, and M = — Tr(a)/2, and 
K = det(a). 3 The covariant Schrodinger equation for particles moving within a thin shell of 
thickness d in 3D is, with both the vector potential V and the electric potential A applied, 1 

iftj^(q,t) = ~G i ^(q ) t) + Q^(q,i) ) (11) 

where Q is the charge of the particle and Dj = Vj — (iQ/h)Aj with Aj being the covariant 
components of the vector potential A. Defining the scalar potential A = —V, we can define 
a gauge covariant derivative for the time variable as D = d t — iQ^oM, and rewrite Eq.(ll) 



as, 1 



ihD Q ij = -^L&DiDrf. (12) 
This equation is evidently gauge invariant with respect of the following gauge transformations: 1 



Aj -> A'j = A, + d jT , A -> A' = A + d tl - ij^^j' = i;e iQ ^ n , (13) 

where 7 is a scalar function. 

Two important facts regarding the wave functions will be needed. 1, the normalization of 
the wave functions remains whatever coordinates are used, and we have with transformation 



of volume element d 3 x =y/Gd 3 q, 1 

J |^(x,t)| 2 d 3 x= J |^(q,t)| 2 V / Gd 3 q=l, (14) 

where 3 ' 4 

G = det(Gy) = g (l - 2Mq 3 + K (q 3 )^ . (15) 

2, an advantage of the curvilinear coordinates is the accessibility of the separability of the 
wave function ?/>(q, t) in (11) or (12) as, 3 ' 1 

t) = Md£l =¥ >(q 3 , t), (16) 

y 1 — 2Mg 3 + K (q 3 ) 2 

and it is guaranteed with suitable choice of gauge for 7 such that A' 3 — 0, 

l(q\q\q 3 ) = - f A 3 (q\ q 2 , q)dq. (17) 
./o 

Combining these two facts, we have two conservations of norm from (14), 

f\x(q\q 2 ,t)\ 2 ^dq 1 dq 2 = l, (18) 



f d 

and / \(p(q 3 ,t)\ 2 dq 3 = 1. 



(19) 



We are now ready to examine the gradient operator V (9) acting on the ip(q,t) and the 
result is, 

M + a 3 K 

V^(q, t) = r^(q, t) + n H - r^x(q\ q\ t)<p(q 3 , t) 

(1 -2Mq 3 + K(q 3 ) 2 ) 1 

+ n ^'^j ^3^(g 3 ,t). (20) 

J 1 - 2Mq 3 + K (q 3 ) 2 

Then taking limit d — > 0, we have, 

V^(q, t) = (r^ + Mn) ^(q, t) + n X (q\ q 2 , t)d q Mq 3 , t), (21) 

which implies that the gradient operator V can be decomposed into two separate parts, one is 
(q l , q 2 ) dependent part (r^d^ + Mn) and another the g 3 -derivative part nd q s, corresponding 
to the decomposition of the Schrodinger equation into two Schrodinger ones determining 
xiq 1 ^ 2 ^) an d y?(? 3 ,£) respectively. 1 Paying attention to the motion on the surface only, 
we have the resultant operator r^d^ + Mn (7). In fact, with proper choice of the confining 



potential V(q 3 ) in the confining procedure, 3 ' 1 we can require that J Q tp*(q 3 , t)d q 3ip(q 3 , t)dq 3 = 
0. So, after performing an integration of operator V in (21) over perpendicular interval [0, d] 
as J Q <p*(q 3 , t)Vip(q 3 , t)dq 3 , only the surface part (r^d^ + Mn) (7) survives. 

The gauge invariance of the momentum operator p = —ihij^d^ + Mn)— QA is assured 
in the presence of the vector potential A with 3D gauge A 3 = being pre-imposed. Under 
2D gauge transformation: A — > A' = A + r M 9 M 7 with 7 = 7(g 1 , q 2 ) and ip — > ip' = e 1 ^ 7 '^, 
we have pip — > p'ip' = e^^pip, 

p'tf,' = {-ih^dp + Mn)-Q(A + r^9 M 7)) ^e iQ7/fi 
= e iQ7/s (-z^(r^a M + Mn)— QA) V 

= e iQ7/s p^. (22) 



So far, we also understand why there is no direct connection between Alb + (M 2 — K) 
and V// + Mn such as in 3D flat space V 2 = V • V. For reaching Alb + (M 2 — K), we have 
to start from the Laplace operator in flat 3D space V 2 = (r M <9 M + nd q s) ■ (r^d^ + nc^) = 
Alb + Md q 3 + 9 2 3 , then resort to the confining procedure. 

geometric momentum IN Dirac's THEORY During 1950's and 1960's, Dirac 11 estab- 
lishes a theory for constrained motion instead follows the routine paradigm of quantization 
hypothesis of kinetic energy T = —h 2 (2fi)ALB on the curved surface. Recalling his famous 
footnote, () we can reasonably infer that if his understanding of canonical quantization is 
self-consistent and indeed insightful, the geometric momentum (5) must be a natural real- 
ization of the momentum in the Dirac's canonical quantization for a system with second-class 
constraints. The second aim of the present study is to illustrate that it is really the case. 

For the constrained motion on the surface 5" (6), Dirac's theory gives for the commutators: 1 

[xi,pj] = ih(5ij - TiiUj), [r,T] = iH—. (23) 

The verification of the first commutator whose tensor form is [r, p] = ih( I — nn) needs 
an identity for the second-rank tensor as I = nn + r M r M whose proof is straightforward. In 
fact, the tensor form of the commutator [r, p] = rp — pr gives: 

[r, p] = -i^[r,(r M ^ + Mn)] = iftr% = ih( I - nn). (24) 
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The second commutator is evident with use of a formula V 2 r = 2Mn, 



[r,T] = —[tAlb] = ^ (-^(^^r)+2(^r)^ 



h 2 



2m 
ih 

— P 

m 



(2Mn+2r%) 



(25) 



GEOMETRIC MOMENTUM DISTRIBUTION OF SPHERICAL HARMONICS The third aim of the 

present study is to give the probability distribution of geometric momentum of the spherical 
harmonics Yi m (6, <p). For our propose, we firstly give GM for particle on the surface of unit 
sphere. 9,16 ' 17 

... „ d sin ip d . . . . 

p x = —inlcos 6 cos ip— -— sm v cos tp), (26) 

dO sm dp 

n d cos ip d . . . . . 

p v = -in{cos0smip— + — — sint/sind, (27) 

oO sm0 dip 

d 

p z = ih(sin 6— + cos 9) . (28) 

These operators satisfy the definition of the vector operator 18 as [Li,pj] = iHEijkPk, we can 
therefore have the operators p x and p y from p z by means of rotation of the axis' rotation. 
Explicitly, rotation 7r/2 around y-axis renders p z to be p x , and — 7r/2 around x-axis renders 
p z to be py, 



p x = exp(-z7rL y /2)p z exp(z7rL y /2), p y = exp(i7rL a ./2)^exp(-i7rL x /2). (29) 

Here we follow the convention that a rotation operation affect a physical system itself. 18 
Hence the eigenvalue problem for operators p x or p y is simultaneously determined once the 
complete solution to p z ^ Vz {ff) = Pzi>p z (9) is known, where on operator p z on the left hand 
side of this equation the carat symbol is used to distinguish the eigenvalue p z on the 
right hand side. The eigenfunctions ip Pz {0) form a complete set once the eigenvalues p z are 
real and continuous, 

^ = ^ tm " P -©- < 30 > 
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They are 5-function normalized, 

^p' z (0, V 2 ) fa. (6, ip) sin 

= — / exp \ i [p , — », m tan - din tan - 
1 /"°° 

= — y exp(?(p' z -p z )z)dz 

= 5(p' z -p z ), (31) 

where the variable transformation In tan 9/2 — y z is used. So, we see explicitly that the 
eigenfunctions il> Pz (6) form a complete set. Next we use it to expand the spherical harmonics 
Yi m (6,ip). Because of the symmetry, the momentum distribution along z-axis depends on 
the angular quantum number I only. The result turns out to be, 

Viipz) = j Y im{0^W Pz (6,<p) sin0d6d<p 

P/(tanh q) 



(32) 



coshg 

where Pi is the Legendre function of order I and the Fourier transform F(p) [f(q)] of a 
function f(q) is defined by, 

fivf-^dq. (33) 

V Z7T 

The first three momentum distribution <pi(p z ) are respectively, 



^ofe) = ^sech(^), (34) 



<Pi(Pz) = 1 -^-p z sech ) , (35) 



^) = ^(3^-l)sech(^). (36) 

The ground state Y 0Q (6,ip) is the minimum uncertainty state for three pairs of (xi,Pi) and 
AxiApi = h/3. In zero angular momentum state Yoo($,ip) that bears no energy either, the 
presence of zero-point the momentum fluctuation (Apj) 2 = H 2 /3 contradicts what classical 
mechanics would indicate. In overall respects, these states bears striking resemblance to 
the probability amplitude of the momentum for one-dimensional simple harmonic oscillator. 
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FIG. 1. Momentum distribution density for Yqq(9,<p) (solid line) and for the ground state of ID 
simple harmonic oscillator (dashed line) are plotted. They are almost identical. 



The momentum distributions in the spherical harmonics offers an experimentally testable 
result for rotational state of spherical molecule such as Cqq. With preparing these molecules 
into ground state of rotation, the probability of the momentum distributions is depicted in 
Fig. 1. 

Conclusions and discussions Bye use of the same confining procedure pioneered by 
Jensen, Koppe 2 and da Costa'^ to give the geometric potential, we find that the momentum 
p = —iftV originally defined in bulk becomes a momentum p = —ih(r fl d fl + Mn) (5) defined 
on the surface, which was previously proposed on completely different ground. Remarkably, 
this momentum is compatible with the Dirac's canonical quantization theory on system with 
second-class constraints. Because r(g x ,g 2 ) (6) in mathematics offers the so-called standard 
parametrization of the 2D surface, the corresponding momentum p (5) should be also prefer- 
able over other forms of momentum such as the generalized momenta (p q i,p q 2) canonically 
conjugated to parameters (g 1 ,^ 2 ). This is another reason it deserves a terminology, geo- 
metric momentum as we called. The distribution amplitudes of the geometric momentum 
of the spherical harmonics are analytically determinable, and experimentally testable for 
rotational state of spherical molecule such as C 60 . 
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